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ABSTRACT 

Thermonuclear bursts from slowly accreting neutron stars (NSs) have proven difficult to detect, yet 
they are potential probes of the thermal properties of the neutron star interior. During the first year of 
a systematic all-sky search for X-ray bursts using the Gamma-ray Burst Monitor (GBM) aboard the 
Fermi Gamma-ray Space Telescope we have detected 15 thermonuclear bursts from the NS low-mass 
X-ray binary 4U 0614+09, when it was accreting at nearly 1% of the Eddington limit. We measured 
an average burst recurrence time of 12+3 d (68% confidence interval) between March 2010 and March 
2011, classified all bursts as normal duration bursts and placed a lower limit on the recurrence time of 
long/intermediate bursts of 62 d (95% confidence level). We discuss how observations of thermonuclear 
bursts in the hard X-ray band compare to pointed soft X-ray observations, and quantify such bandpass 
effects on measurements of burst radiated energy and duration. We put our results for 4U 0614+09 in 
the context of other bursters and briefly discuss the constraints on ignition models. Interestingly, we 
find that the burst energies in 4U 0614+09 are on average between those of normal duration bursts 
and those measured in long/intermediate bursts. Such a continuous distribution in burst energy 
provides a new observational link between normal and long/intermediate bursts. We suggest that 
the apparent bimodal distribution that defined normal and long/intermediate duration bursts during 
the last decade could be due to an observational bias towards detecting only the longest and most 
energetic bursts from slowly accreting NSs. 

Subject headings: X-rays: bursts — binaries: close — X-rays: individual (4U 0614+09) — stars: 
neutron — X-rays: binaries — accretion, accretion disks 



1. INTRODUCTION 

Matter accreted onto neutron stars (NSs) in low-mass 
X-ray binaries (LMXBs) piles up and settles into the NS 
envelope. After a period accumulating fuel, thermonu- 
clear reactions become unstable and the accreted layer 
ignites, producing a thermonuclear flash (also known as 
type I X-ray burst, "burst " here after; Grindla v et al.l 
Il976t iMaraschi fc Cavalierd 119771 ). During the ther- 
monuclear runaway all or part of the accreted hydro- 
gen (H) and/or helium (He) are converted into heavier 
elements so that fresh fuel needs to accumulate before 
the next burst occurs, giving rise to a cyclic behavior. 
This has been observed in many thermonuclear burst 
sources ("bursters") that alternate short bursts (typi- 
cally seconds- to minutes-long) with extended periods of 
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persisten t accretion-powered emission (typically hours- 
days; see lLewin et al.fi"9 93; Stro hmayer fc Bildstenll200"6l 
for reviews). The persistent luminosity, L pers , often ex- 
pressed as a fraction of the Eddington limit, is thought 
to trace the mass accretion rate per unit surface area, 
rh. When the rh history is known, the time spanned be- 
tween two successive bursts (known as burst wait time, 
^wait) determines how much mass needs to be piled up 
until the ignition conditions are met at the base of the 
accreted layer. When data gaps prevent the detection 
of successive bursts a more convenient parameter is the 
wait time averaged over the available observations, i.e., 
the ratio between total accumulated exposure time and 
number of bursts detected, which we refer to as burst 
recurrence time, i rcc . 

The vast majority of thermonuclear bursts observed 
during more than three decades, referred to as normal 
bursts, last between a few seconds and a few minutes 
and r adiate between 10 38 e rg and a few times 10 39 erg 
(e.g., Gallowa y et al.1 12008? ) . In the last decade two 
new classes of bursts have emerged: long/intermediate 
bursts (long duration bursts or simply "long bursts" ; e.g., 
lin't Zand et al.ll2005|) , which last tens of minutes and ra- 
diate 10 40 ~10 41 erg; and superbursts, which last for hours 
up to more than a day and radiate energies of the or- 
der of 10 42 erg (ICornelisse et al.l [20001: IKuulkersI [200l 
iChenevez et~aLl I2008D . ~As" more long bursts have been 
discovered, the question was raised of whether they form 
a distinct population of bursts or they constitute a rare, 
extremely long and energetic case of normal bursts. 

Long bursts take place preferentially when rh is low, 
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near or below 1% of the Eddington rate, thus their prop- 
erties depend on the nuclear burning processes at work at 
the lowest accretion rates, where sedimentati on of heavy 
elements may play a role (jPeng et al.l 120071 ) . Further- 
more, ignition of long bursts takes place deep in the en- 
velope after a thick layer of fuel is built up. Hence their 
properties are also sensitive to the heat flowing from the 
crust, providing an independent meth od to study the 
therm al properties of the NS interior (|Cumming et al] 
2006). We refer hereafter to bursts occurring when m 
is close to or lower than 1% of the Eddington limit as 
low-m bursts, with the understanding that these include 
both normal and long duration bursts. 

4U 0614+09 harbors a slowly and persistently ac- 
cretin g NS (active since its discovery, Gi acconi et al.1 
ll972Q and belongs to the sub-class of low lu- 
minosity NS-LMXBs (the so-called atoll sources, 
iHasinger &: van der Klisl Tl989f ). The system is likely 
an ultraco mpact binary, based on the optical-to-X-ray 
flux ratio (Ivan Paradiis fc McClintockl 11994 Uuett et al.l 
20011). tentative orbit al periodicities ( O'Brien et al.ll2005F 



Shah baz et all I2008D and its low but persistent X-ray 
luminosity (jin't Zand et all 120071 ). However, the ex- 
act nature and composition of the (likely degenerate) 
donor star has been a matter of debate (although strong 
constraints ha ve been placed from X-ray and optical 
spectroscopy; Uuett et all 120011 : iNelemans et al.1 120041 : 
iWerner et al.ll2006l ). 

Although 4U 0614+09 is a known burster (jSwank et al.1 
Il978t iBrandt et al.l 119921 ) and a frequently observed 
source, a relatively small number of bursts have been 
recorded from 4U 0614+09 since its discovery four 
decades ago. Only two bursts were detected with the 
Rossi X-ray Timing Explorer (RXTE) despite numerous 
observations over the course of 16 yr (which amount to 
more than 2 Msec of on-source exposure time). A com- 
pilation of these and 31 other bursts det ected with sev- 
eral X - ray o bservatories was presented by Kuulkcr s et al.1 
(|2010l IKIOI hereafter). They found two long duration 
bursts during 2001-2002 with integrated radiated ener- 
gies E'b >4xl0 40 erg, when the persistent soft X-ray lu- 
minosity was steady, and an increased normal burst rate 
after that period with the corresponding t rec in the 9- 
14 d range. 

The scarcity of bursts detected from 4U 0614+09 or 
similar systems is not surprising. Low-m bursts are ex- 
pected to recur on timescales much longer than the ob- 
serving windows of pointed X-ray telescopes, which typ- 
ically observe LMXBs with little or no interruptions for 
up to about a day. Hence detection with pointed in- 
struments is rare and both i wa ;t and t rcc have remained 
poorly constrained. Thanks to its very large instanta- 
neous field of view (~10 sr) and its sensitivity to X-ray 
photon energies down to 8 keV, th e Gamma-ray Burst 
Monitor (GBM; iMeegan et"aTll2009D on board the Fermi 
Gamma-ray Space Telescope is able to detect rare, bright 
and short-lived X-ray bursts. This is of particular rel- 
evance to the study of thermonuclear bursts from ac- 
creting NSs. A dedicated search using wide-field X-ray 
monitors provides the most efficient way to measure ac- 
curately the recurrence time of low-m bursts. With this 
goal, we started on 2010 March 12 a systematic GBM 
all-sky search for X-ray bursts, as described in Section [2j 
We report 15 thermonuclear bursts from 4U 0614+09 



collected during the first year of our campaign, between 
2010 March 12 and 2011 March 12. We present the burst 
properties and accretion rate history in Section[3] In Sec- 
tion|4]we discuss these results, quantify bandpass effects 
on the observed burst properties (Section I4.1j) and com- 
pare 4U 0614+09 to other low-m bursters (Section 14. 2|) . 
Section [5] summarizes our main findings and conclusions. 

2. DATA ANALYSIS 

2.1. GBM all-sky search for X-ray bursts 

GBM is an all sky monitor whose primary objective is 
to extend the energy range over which gamma-ray bursts 
(GRBs) a re observed in the L arge Area Telescope (LAT) 
on Fermi (jMeegan et al.ll2009l ) . The twelve sodium iodide 
(Nal) scintillation detectors cover an energy range from 8 
keV to 1 MeV. Due to the spacecraft's 50 degree rocking 
angle, exposure of the sky is essentially uniform through- 
out the current search period between 2010 March 12 
and 2011 March 12. The X-ray burst search uses the 
first 3 channels (8-50 keV) of CTIME data from the Nal 
detectors. Periods where the spacecraft is performing 
a rapid maneuver (>2.75xl0 -3 rad/s), SAA passages, 
GRBs and contaminating features, such as solar events, 
that prevent a good background model fit are removed 
from the data. An empirical background model which 
includes Earth occultations of bright sources is fit to the 
data. Short-lived (~10-1000 s) X-ray flares are selected 
by visually inspecting the CTIME Channel 1 data (12-25 
keV) along with the fitted background model. Using the 
first three channels, a preliminary localization is obtained 
using the cosine differential response of the detectors (see 
Sec. 12. 2p . If the \ 2 °f the localization is less than 450, 
the localization is more than 10° or 3cr away from the 
Sun and the position lies above and more than la from 
the Earth's limb (a 66° radius circle around the geocen- 
tcr), the selected feature is considered an X-ray burst 
candidate and set aside for further analysis. We define 
the observing duty cycle as the ratio between the time 
when the source was visible to GBM (excluding times of 
Earth occultation, rapid slews and SAA passages) and 
the total time spanned by the current search (one year). 
In the case of 4U 0614+09, the duty cycle was 49.3% and 
the total effective exposure time 180 d. 

In order to optimize the signal-to-noise ratio (S/N) 
in the X-ray bursts from 4U 0614+09 and measure 
the burst timescales and peak count rate, we extracted 
background-corrected 12-25 keV light curves adding the 
three brightest detectors, using daily CTIME data re- 
binned to Is time resolution. We defined the burst rise 
time £ r i sc as the time for the resulting intensity to increase 
from 25% to 90% of its peak value. We interpolated lin- 
early between data points and we estimate an uncertainty 
of 0.5s (half the time resolution used) in f r ; se . We used 
CSPEC light curves in the 8-15 keV band keeping the 
original 4s time resolution in order to measure the du- 
ration of the bursts, defining the burst end as the time 
when the intensity drops below 10% of the peak. 

2.2. Localization 

The 12 Nal detectors are thin disks with different ori- 
entations. Using the known angular response of the de- 
tectors, the relative count rates registered during a burst 
in the 12 detectors provide a solution for the likely ar- 
rival direction of the burst. The measured count rates 
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TABLE 1 

GBM X-RAY BURSTS FROM 4U 0614+09. 



ID 




Peak time 




Dctcctors a 


Peak rate b 


rt. A. 




Error 






(UTC) 






n ( ) 


(c/s) 


( ) 


( ) 


( 1 


Bl 


2010- 


-03-28 00: 


:16: 


:13 


9(22) 0(42) 1(44) 10(45) 2(59) 


381+37 


98.4 


4.6 


3.4 


B2 


2010- 


-04-29 03: 


:46: 


:49 


1(22) 2(26) 0(45) 5(59) 


617±41 


92.6 


8.5 


6.0 


B3 


2010- 


-05-18 10: 


:56: 


:59 


4(29) 5(32) 3(50) 


360+37 


95.0 


4.4 


9.3 


B4 


2010- 


-05-25 07: 


:35: 


:59 


1(14) 0(19) 3(46) 5(57) 2(59) 


595±40 


104.6 


1.4 


6.8 


Bo 


2010- 


H7 ncr nc: 
-Ui-UO uo. 


.Uo. 


:03 


5(16) 1(43) 2(51) 3(55) U(60) 


593+41 


87.8 


6.6 


5.5 


Bfi 


2010- 


-07-29 07: 


:12: 


:47 


1 (26) 2(29) 5f39i 0(^0) 


496+39 


99.0 


10.5 


8.2 


B7 


2010- 


-08-09 09: 


:30: 


:02 


4(17) 3(36) 5(46) 


740+43 


87.9 


10.1 


2.9 


B8 


2010- 


-09-10 18: 


:20: 


:15 


3(14) 4(33) 5(60) 


489+35 


96.7 


15.2 


4.5 


B9 


2010- 


-10-12 03: 


:54: 


:43 


10(37) 2(37) 9(46) 1(47) 


311+34 


94.9 


13.7 


6.2 


BIO 


2010- 


-10-21 09: 


33: 


:42 


7(23) 6(32) 3(42) 8(46) 4(60) 


596+41 


94.3 


9.6 


3.5 


Bll 


2010- 


-11-28 03: 


:50 


:45 


9(27) 10(30) 11(41) 


668+40 


97.0 


12.0 


4.1 


B12 


2011- 


-01-11 20: 


:29: 


:51 


9(18) 11(44) 10(45) 6(51) 7(54) 


490+39 


87.0 


6.8 


3.6 


B13 


2011- 


-01-14 15: 


:26: 


:50 


9(27) 11(36) 10(40) 7(59) 


302+36 


97.7 


18.4 


9.7 


B14 


2011- 


-01-18 22: 


■AT 


:38 


7(12) 6(17) 9(50) 0(57) 8(57) 11(60) 


627+41 


98.5 


13.2 


10.4 


B15 


2011- 


-01-24 11: 


:00: 


:23 


9(22) 10(31) 11(46) 


716+42 


89.8 


9.4 


7.8 



a List of detectors pointing within 60° of 4U 0614+09 at the time of the burst, used for spectral analysis (Sec. 1273)1 , The angular distance 
from each detector to 4U 0614+09 is given between parenthesis, in degrees. Only detectors with pointing offsets <50° were used to extract 
the light curves. 

b Net 12-25 keV count rate from the three brightest detectors. 



are compared with model rates for arrival directions on 
a 1° resolution grid measured in the spacecraft coordi- 
nate frame. A ^-minimization algorithm finds the best 
match for the observed rates to 41168 positions on the 
sky, with an uncertainty that depends on how steeply 
\ 2 increases away from the minimum. In order to make 
the assessment of the goodness-of-fit independent of the 
intensity of the burst, the observed rates are first normal- 
ized to a fiducial "average" burst before calculating x 2 - 
This reduces the % 2 values for very bright bursts with 
low statistical uncertainties where the systematic errors 
dominate, and allows us to use one \ 2 value to reject 
poor fits to the data over the range of burst intensities 
detected by GBM. 

Because the angular response is energy-dependent, a 
source spectrum is folded through the response when cal- 
culating the model rates. The model rates in the search 
for X-ray bursts are calculated at energies below 50 keV 
and compared to count rates measured in the correspond- 
ing CTIME channels. To maximize the signal in the 
search for these events, the counts in CTIME channels 
through 2 are used, and the model rates are evaluated 
between 8 and 50 keV. The uncertainties in the detector 
responses below 10 keV and the noisier background at 
the lowest energies may introduce larger systematic er- 
rors than incurred above 10 keV, so the localization is 
also performed in the 10-50 keV energy band. By de- 
fault the search for X-ray bursts uses a power-law model 
rate function with an index of 2, which represents a va- 
riety of phenomena but may not be the best model for 
type I X-ray bursts. Model rates for a softer power-law 
function with index 3 were therefore used at the localiza- 
tion stage once the X-ray burst candidates were isolated 
in the overall data set. Model rates with a blackbody 
spectrum that was later shown to be a good fit to the 
X-ray burst candidates were subsequently generated and 
used to see whether a good model fit produced a better 
localization. 

As expected, the statistical error on the localization 
was higher by about 1° for a given model when the 
rates below 10 keV were excluded from the localization 



of the 15 events associated with 4U 0614+09. The er- 
ror on the localization, inferred from the distance to the 
known source position, was comparable for localizations 
performed using both energy ranges. The ratio between 
the offset from the source and the statistical uncertainty 
suggests that the systematic error on the localizations 
using the 8-50 keV energy range may be larger, with 
41% (82%) contained within the l(2)<r uncertainty region 
compared to 54% (87%) for the localizations performed 
in the 10-50 keV energy range. The centroids (weighting 
each position by its statistical uncertainty) of the 15 com- 
bined locations are between 0.5°-l.l° from 4U 0614+09 
for the 10-50 keV range localizations depending on the 
model, with the more representative softer power-law and 
blackbody spectral models yielding centroids at 0.5 and 
0.6° from the source. The blackbody and softer power- 
law models produced localizations that were on average 
closer to the source by about 0.7° compared to the power- 
law model with index of -2. The 8-50 keV localizations 
produce a centroid that lies between 1.6 and 1.9° from 
the source. Within the relatively low (15) number of 
events available, our results suggest that systematic er- 
rors may be larger when the full energy range is used to 
localize events. Whilst our analysis suggests that using 
the 10-50 keV energy range with either a soft power-law 
or black-mody model produces the best results when lo- 
calizing type I X-ray bursts using GBM data, the events 
are weak enough that the uncertainties are dominated 
by statistics, and these systematic effects arc small by 
comparison. 

2.3. Spectral analysis 

We performed time-averaged and time-resolved spec- 
troscopy of the 15 confirmed bursts from 4U 0614+09. 
We used the continuous CSPEC data type in our analy- 
sis, with the best availa ble spectral resolutio n for GBM 
and 4 s time resolution (Mcega n et al.ll2009l ). The data 
were analyzed with the RMFIT (3.3rc8) software pack- 
age which has been developed for the GBM data analysis. 
The time- varying background was fitted with polynomial 
functions of the third or fourth order. We generated De- 
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tector Response Matrices using GBMRSP VI. 9. In our 
analysis we used all the Nal detectors with zenith an- 
gles to the source smaller than 60°, and which view was 
not obstructed by the solar panels or other parts of the 
spacecraft. 

We have fitted the time-integrated spectra with a 
power law and a blackbody function, using the Castor 
statistic (C-stat m. and in all bursts a blackbody pro- 
vided the best fit to the data. The difference in C-stat 
between the two models ranges from several tens to hun- 
dreds, illustrating the improvement of a thermal model 
over a power law. To quantify this and estimate the 
significance of the blackbody fit improvement over the 
power law fit, we simulated 10 4 spectra from the best- 
fit model parameters. We found that for a typical value 
(~40) of the difference in C-stat between thermal and 
power law model fits, the fit improvement is significant at 
the >99.99% (3.9er) confidence level. Since these events 
are spectrally very soft for the GBM energy band, we 
used the 8-30 keV range in our analysis, but we note that 
using the full Nal energy range gives consistent results. 
We also checked spectral fits from individual detectors 
and they all give consistent results. Here we use the re- 
sults obtained by combining the appropriate detectors to 
maximize the S/N. 

2.4. Persistent emission and mass accretion rate 

In order to characterize the soft and hard persistent 
(accretion-powered) X-ray emission from 4U 0614+09 
and estimate m, we used 1.5-12 ke V light curves col- 
lected by the All-Sky Monitor (ASM; iLevine et ai1ll996D 
on board the Rossi X-ray Timing Explorer as well as 15- 
50 ke V light curves from Sw ift's Burst Alert Telescope 
(BAT; iBarthelmv et al.l feOOSl , using data fro m the hard 
X-ray transient monitor ([Krimm et al.l 1200 6;) . We used 
ASM and BAT daily intensity measurements taken be- 
tween 2010 March 12 and 2011 March 12 including all 
available (2-cr) detections, which amount to 81% and 
76% of the whole one-year period for ASM and BAT, 
respectively. 

As we are interested for the purpose of this work in 
long-term (weeks to months) changes in m as well as its 
mean value over the search period, we applied a running 
average with a 12-d window to the resulting light curves. 
At the expense of disregarding variability on timescales 
of less than ~2 weeks, this method increases the S/N and 
removes short data gaps due to, e.g., high background in 
ASM or BAT. The window length was choosen to pre- 
vent data gaps from affecting the running average, while 
keeping a relatively short window. We then converted the 
ASM and BAT intensity into 2-20 keV and 20-100 keV 
luminosity, respectively, using WebPimms3 and a dis- 
tance to 4U 0614+09 of 3.2 kpc ([Kiel ). We assumed 
that the X-ray spe ctrum of 4U 0614+09, averaged over 
its spectral states (|Linaresl 12009). is well represented by 
a power law model with photon index 2. The resulting 
2-100 keV luminosity gives an estimate of the bolometric 
persistent luminosity, £ pcrs , which we convert into m as- 
suming that all the gravitational energy is radiated away 

11 A modification of Cash's C-statistic used to fit spectra with 
low number of counts, which tends to y 2 for large number of counts; 
http: / /heasarc. gsfc.nasa.gov/docs/xanadu/xspec/wstat.ps 

lz Portable Interactive Multi-Mission Simulator v. 4.5 at 
http: / /heasarc. gsfc.nasa.gov/Tools/w3pimms. html 



isotropically at the surface of a 10 km radius, 1.4 M Q 
mass NS (this and other assumptions may result into 
un certainties in m o f a factor ^2; see, e.g., discussion 
in iCoriat et al.1 l2012t ) . When normalizing L pcrs by the 
Eddington luminosity (^Edd) we use a fiducial value of 
2.5xl0 38 erg s _1 , yet we note that the exact value de- 
pends on the composition of the accreted material. Using 
3.8xl0 38 erg s _1 , a value more typical of H-poor mate- 
rial, does not change our results and conclusions. 
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55300 55400 55500 55600 
Time (MJD) 

Fig. 1. — Bottom panel: Evolution of mass accretion rate per 
unit surface (left axis) derived from the 2—100 keV persistent lu- 
minosity (rig ht a xis) after applying a 12-d window running aver- 
age (see Sec. 12.41 for details). Red arrows along bottom axis show 
the times of the 15 thermonuclear bursts from 4U 0614+09 de- 
tected with GBM during the 1 yr period. Month start dates are 
marked along the top axis, from March 2010 to March 2011. Top 
& middle panels: Hard (BAT) and soft (ASM) X-ray light curves 
of 4U 0614+09. Daily-averaged intensity measurements (left axes) 
are converted into luminosity (right axes). 



3. RESULTS 

We found a total of 15 X-ray bursts from 4U 0614+09 
during the first year of our Fermi-GUM X-ray burst 
search (Sec. l2.ip . Only 33 bursts had been observed from 
4U 0614+09 since it was discovered, using observations 
taken over 32 yr by 8 different X-ray observatories (be- 
tween 1975 and 2007; see|K10|). Therefore, the 15 GBM 
bursts presented herein constitute a substantial improve- 
ment in burst detection rate. Table [T] shows all GBM 
bursts labeled in chronological order B1-B15, their peak 
times, count rates and locations. All bursts were clearly 
detected in the 12-25 keV band with all GBM Nal de- 
tectors pointing within 60° of 4U 0614+09, and none of 
them was detected above 50 keV, as expected for ther- 
monuclear bursts given their thermal spectrum. Figure[T] 
shows the burst times as well as the soft (ASM) and hard 
(BAT) X-ray light curves of 4U 0614+09. From these, 
we obtained the 2-100 keV persistent luminosity and in- 
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TABLE 2 

Properties of the GBM type I X-ray bursts from 
4U 0614+09. 



ID 


kT poak 


Lpcak 


E 


Rise 


Duration 


At 




(koV) 


(10 3S crg/s) 


(10 39 erg) 


oo 


(s) 


(d) 


Bl 


2.7±0.1 


1.79±0.07 


2.07±0.01 


1.6 


28.6 




B2 


3.0±0.1 


1.65±0.05 


1.77±0.02 


2.2 


20.3 


32.1 


B3 


2.9±0.2 


0.89±0.07 


0.72±0.03 


0.9 


29.4 


19.3 


B4 


2.9±0.1 


1.75±0.05 


1.06±0.03 


1.2 


11.1 


26.2 


B5 


2.6±0.1 


1.97±0.07 


2.33±0.12 


1.5 


26.0 


40.9 


B6 


3.1±0.1 


1.28±0.06 


0.82±0.02 


1.1 


16.4 


24.4 


B7 


3.2±0.1 


2.72±0.06 


5.77±0.03 


4.5 


46.5 


35.2 


B8 


3.0±0.1 


1.91±0.05 


4.19±0.03 


1.8 


44.0 


32.4 


B9 


2.7±0.1 


1.38±0.06 


2.23±0.02 


1.3 


34.1 


31.4 


BIO 


3.0±0.1 


2.36±0.09 


2.15±0.06 


1.5 


20.7 


9.2 


Bll 


3.2±0.1 


2.17±0.07 


2.84±0.03 


2.8 


30.5 


37.8 


B12 


3.0±0.1 


1.89±0.08 


2.58±0.05 


0.9 


49.0 


44.7 


B13^ 


3.3±0.2 


0.74±0.05 


0.43±0.02 


1.3 


10.1 


2.8 


BlT 


2.9±0.1 


2.05±0.08 


1.49±0.07 


2.9 


14.0 


4.3 


B15 


3.2±0.1 


2.30±0.08 


1.82±0.08 


3.1 


15.0 


5.5 



a Burst peak luminosity L pea k and radiated energy E at 3.2 kpc. 
Rise time, duration and E as measured by GBM, not corrected for 
bandpass. At is the time in days since the previous burst from 
4U 0614+09 detected by GBM (not corrected by on-source expo- 
sure time and therefore corresponds to an upper limit on the burst 
wait time, t W ait; Sec. [TJ. 

b Due to the low number of counts per spectrum, blackbody tem- 
perature and radius are strongly correlated in B13. If we fix the 
radius to a value of 5 km, consistent with all other bursts, the peak 
temperature becomes kTp,,^ ~2.5 keV. 
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Fig. 2.— GBM localizations of the 15 bursts from 4U 0614+09 
presented i n th is work (gray dashed circles show la error circles; 
see Section E2] and Table Q) . The average position from all GBM 
bursts (blac k open circle), the k nown position of 4U 0614+09 (black 
filled circle; Migliari ct al. 2006) and the closest bright X-ray source 
(Crab, black diamond) are also shown. 



ferred mass accretion rate per unit area m (Sec. 12 .4[) . 
which remained between 0.6% and 1.3% of the Edding- 
ton limit (LEdd and rn-Edd, respectively; Sec. 12 .4[) during 
the whole period studied in this work, with an average 
value < Lp era >= 2.3xl0 36 erg s _1 (m/m Edd = 0.9%). 
All the RXTE observations taken between March 2010 
and March 2011 showed 4U 0614+09 in the intermedi- 



ate state, where it remains most of the time (the accre- 
tion state inter mediate betwee n spectrally hard and soft 
states; see, e.g., Linares 2009). 

As shown in Figure [2j the burst locations are consis- 
tent with the known position of 4U 0614+09. For about 
90% of the bursts the 2a confidence region around the 
best location includes 4U 0614+09 (see full details in 
Section |2~2|) . All bursts are more than 10° and 70° away 
from the center of the Sun and Earth, respectively, as 
seen from the Fermi spacecraft (Sec. 12.11) . The closest 
known bright X-ray source (Crab) and burster (4U 0513- 
40) are about 16° and 51° away, respectively, and there 
is no known LMXB within 15° of 4U 0614+09. This, to- 
gether with the thermonuclear nature of the bursts (see 
below), allows us to identify all bursts B1-B15 as coming 
from 4U 0614+09. 
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Fig. 3. — Time-resolved spectroscopy of the GBM X-ray bursts 
from 4U 0614+09, showing (from top to bottom): net count rate, 
bolometric flux, blackbody temperature and blackbody radius. 
The temperature decay is clearly visible, identifying all bursts as 
thermonuclear. 



We first performed power law fits to the time-averaged 
8-30 keV spectrum of all bursts and obtained photon in- 
dices in the 3-4.4 range. Such high values of the photon 
index reveal a relatively steep or "soft" spectrum, in- 
dicative of thermal emission. Indeed, a blackbody model 
gives a better fit (Sec. I2.3|) to the average spectrum of all 
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Fig. 4. — Bottom: Measurements of burst recurrence time t rGC in 4U 0614+09, between 1996 and 2012. The red filled rectangle shows 
our measurement of t rec from the first year of the GBM all-sky X-ray burst search. The red vertical error bar marks the shortest wait time, 
2.8 d , and the crossed red rectangle shows the t Tec measured in January 2011. Gray dashed rectangles display the limits on t rC c found by 
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bursts, yielding temperatures in the 2.5-2.8 keV range, 
fully consistent with a thermonuclear origin. The results 
of the time-resolved spectroscopy (Sec. 12.31) are shown in 
Figure |3l and confirm the thermonuclear nature of the 
bursts. We detect in all cases a significant drop in burst 
temperature along the decay, also referred to as "cooling 
tail" . This is a defining property of type I X-ray bursts 
and a sufficient condition to identify bursts as thermonu- 
clear. Burst peak temperatures (ATp^t) are in the range 
2.6-3.3 keV (Table©, and decay to ~2 keV in the latest 
detected stages of the bursts. 

We use thro ughout this work the 3.2 kpc distance 
found by IKlOl from the brightest photospheric radius 
expansion (PRE) burst in their sample. The apparent 
blackbody radii that we find (Figure |3] between ~3 km 
and ~10 km during most phases of the bursts) are con- 
sistent with previously reported values for 4U 0614+09 
(|K10l ) and compatible with the expected NS emitting 
area (note that no gravitational redshift or spectral cor- 
rections have bee n applied to these values; see, e.g., 
iLewin et al.l [1993) . All the bursts we detect feature 
single-peaked lightcurves, also when inspecting 100-200 s 
around the peak with 1-4 s time resolution, and we 
find no signs of PRE in the GBM time resolved spec- 
troscopy of any of the bursts (but see Sec. I4.1|) . The 
burst peak luminosities (L pca k) that we find range from 
0.7xl0 38 erg s" 1 to 2.7xl0 38 erg s" 1 (Tabled]). The rise 
times we measure in the GBM band (t r i se ; Sec. 12.11 for 
definition) are between 1 s and 5 s. By integrating the 



bolometric luminosity over the fraction of the bursts de- 
tected by GBM, we find burst energies in the range [0.4- 
5.8] xlO 39 erg. The burst duration in the GBM band was 
between 10 s and 50 s. Due to the combined effects of the 
GBM bandpass and background, these are lower limits 
to the total radiated burst energy and duration (see Sec- 
tion 14.11 for further discussion of bandpass corrections) . 

From the 15 bursts detected and the ~50% observing 
duty cycle (Sec. 12. 1|) . we find a burst recurrence time 
between March 2010 and March 2011 of t rec =12±3 d 
(la confidence interval from Poisson statistics). Notably, 
bursts B12 and B13 were only 2.8 d apart, the closest 
burst pair ever detected from 4 U 061 4+09 (shortest t W ait 
known to date was one week; IKlOl ). Four bursts were 
detected in January 2011 corresponding to a recurrence 
time of ~4 d, shorter than the recurrence time during 
the preceding 10 months (i rcc was ~13 d between 2010 
March 12 and December 31). This can be seen in Fig- 
ure [TJ which shows the times of the bursts detected by 
GBM together with the m evolution. Assuming a Poisson 
distribution, the probability of detecting at least 4 bursts 
during one month (effectively two weeks) for the aver- 
age burst detection rate 1.3 month -1 (effectively 1.3 per 
2 week-period) is 3.8%. This suggests that the intrin- 
sic thermonuclear burst rate from 4U 0614+09 increased 
during January 2011. The source state during January 
2011 was similar to the rest of the search period. The 
inferred m reached a minimum near January 2011 (at 
w/rh-Edd — 0.6%). 
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Fig. 5. — Spectral evolution during a long PRE burst from 4U 0614+09 detected in 2001 (open gray circles; from lKlOl . burst FR 51944). 
Panels show, from top to bottom: bolometric flux, blackbody temperature and radius (for a distance of 3.2 kpc and without color or 
gravitational redshift corrections). We compare this to the spectral evolution measured by GBM during burst B7, detected in 2010 (filled 
red squares; see Tables [JJ & |2j . Left panels display the entire (about 8000 s long) long burst time range. Right panels zoom into the 
touchdown phase of the burst, the only part hot enough to be detected in the hard X-ray band (i.e., by GBM and FREGATE). 



We do not detect long/intermediate bursts from 
4U 0614+09 during the 1 yr search period, and we can 
place a 95% lower limit (assuming a Poisson distribu- 
tion, iGehrelsl [1986) on the long burst recurrence time 
of 62 d. Our measurement of t Iec is consistent with 
previous estimates, which placed t rsr between one week 
and t hree months (jin't Za nd et al.ll2007t iKuulkcr s et al.l 
120101) . Figure U compares t Tec from 4U 0614+09 as 
measured by several instruments during the last sixteen 
years, between 1996 and 2012. In only one year, the 
GBM results presented herein have provided the most 
accurate measurement of t Tec and revealed the shortest 
iwait in 4U 0614+09. Furthermore, we are beginning to 
probe changes in t rcc on timescales of months, as illus- 
trated by the likely increase in burst rate during January 
2011. 

4. DISCUSSION 

We have presented the detection and properties of 15 
thermonuclear bursts from 4U 0614+09 observed be- 
tween March 2010 and March 2011 within our GBM 
all-sky search for X-ray bursts. Using only one year of 
GBM data we increased by 40% the sample of bursts 
from 4U 0614+09 collected over more than three decades 
(K10). These results illustrate the benefits of an X-ray 
monitor with instantaneous coverage of most (~75%) of 
the sky and yield a robust measurement of the recurrence 
time of low-m bursts in 4U 0614+09: i rec =12±3 d. 

Figure H shows the full 1996-2012 ASM light curve 
(top panel; monthly-averaged intensity with error bars 
showing its standard deviation during each month) to- 
gether with i rec in 4U 0614+09. The average burst rate 
has been steady over 16 yr, with two marginal ch anges : 
a decrease in burst rate in 2001-2002 reported by IKIOI 
associated with lower variability in the ASM light curve, 
and the increased burst rate in January 2011 reported 



herein. The two only known long/intermediate bursts 
from 4U 0614+09 happened in the 2001-2002 period, 
when the variability on long (days-months) timescales 
was at a minimum (Fig. |4]). This "calm" period has not 
repeated ever since, which could explain why we do not 
detect long/intermediate bursts during 2010-2011. 

4.1. Bandpass corrections: total burst energy and 
duration 

The energy spectrum of thermonuclear bursts is for 
most purposes well described by a simple absorbed black- 
body model. At the peak (near the flux maximum) 
bursts can show a wide variety of light curve shapes 
and spectral evolution. In particular, the peak of photo- 
spheric radius expansion (PRE) bursts features a dras- 
tic increase of the blackbody radius Rbb and an asso- 
ciated drop of the blackbody temperature fcTbb, due to 
the effect of radiat ion forces on the NS photosphere (e.g., 
iLewin et al.lll993l and references therein). Subsequently 
the photosphere contracts and heats up until it reaches 
the NS surface again, at "touchdown". Most low-m 
bursts are PRE events (of either normal or long dura- 
tion class) as the high ignition depths power energetic 
and luminous events that can overcome the gravitational 
force at the photosphere. 

Soft X-ray sensitivity is necessary in order to observe 
the full radius expansion phase of PRE bursts, due to 
the simultaneous drop in temperature which shifts the 
spectrum to lower energies. Conversely, the lack of spec- 
troscopic evidence for PRE when observing a given burst 
in the hard X-ray band does not rule out the PRE na- 
ture of the burst. This is shown in Figure [5j where 
we compare the spectral evolution of a long PRE burst 
from 4U 0614+09 observed in 2001 with broadband spec- 
tral coverage (wi th th e PCA, ASM and FREGATE in- 
struments; from IKIOD to one of the GBM bursts pre- 



TABLE 3 

Bandpass corrections to burst energy and duration. 



Burst 



GBM 
(8-1000 keV) 
Eobs (erg) tobs (s) 



3.6xl0 39 
3.4xl0 40 



FREGATE 
(6-400 keV) 
-Bobs (erg) fobs 



6.0xl0 39 
3.9xl0 40 



PCA 
(2-60 keV) 
E b (erg) t h (s) 



FREGATE/GBM 



PCA/GBM 



-^obs /-^obs ^obs /^obs -^b/-^obs /^obs 



Normal (FR 52951) 
Long (FR 51944) 



30 
140 



15 

308 



? 



4.4 xlO 40 



? 

7700 



1.7 
1.1 



1.5 
2.2 



>1.7 
1.3 



>1.5 
55 



a Burst energy and duration for a lo ng an d normal burst measured with three instruments in different ener gy bands, as indicated. PCA 
and FREGATE data are taken from iKlOt GBM values are obtained from simulated light curves (Sec. |4~TT |. PCA values are considered 
bolometric measurements of the burst energy and duration. Bandpass correction factors are given as ratios between the values measured 
in: FREGATE/GBM and PCA/GBM. Based on the last two col umn s we estimate the bolometric corrections on burst energy and duration 
as Eb/^obs = [l-3-1.7] and i bsAobs = [1-5-55], respectively (Sec. |4~TT |. 



sented herein. In all cases except B3 and B13 the fluxes, 
temperatures and blackbody radii measured during our 
GBM bursts are consistent with those measured near the 
"touchdown" phase of a PRE burst (Figure [5] shows one 
example, B7). Thus most thermonuclear bursts detected 
with GBM (and perhaps with any hard X-ray wide field 
detector) are in fact consistent with PRE bursts in the 
touchdown phase. However, as our GBM X-ray burst 
search includes faint and untriggered events we are also 
sensitive to fainter, sub-Eddington bursts, as witnessed 
by the detection of B3 and B13 (Table [2] & Figure [3]). 
All the GBM bursts presented herein are much shorter 
than the segment of the long burst detected by FRE- 
GATE. This, together with the simulations described be- 
low, indicates that all bursts B1-B15 are normal duration 
bursts, even after bandpass corrections are applied. 

Due to its spectral evolution, a given burst can 
produce very different light curves when observed 
with instruments with different bandpasses (see, e.g., 
IChelovekov eTail [2006h . This is of particular relevance 
when studying thermonuclear bursts with a hard X-ray 
detector, as many of the burst photons are emitted out- 
side the energy range to which the instrument is sensi- 
tive to. Whenever a burst is clearly detected, its thermal 
spectrum is usually well constrained using the hard X-ray 
band and the bolometric flux can therefore be measured, 
under the assumption that color and effective tempera- 
tures are equal. 

However, the actual burst detection relies on having 
sufficient burst photons emitted in the instrument band 
to rise significantly above the background rate. In prac- 
tice, this implies that instrum ents like GBM (sens itive to 
photon energies above 8 keV; lMeegan et al.ll2009l) FRE- 



GATE (sensitive above 6 keV; lAtteia et al.l 120031) ■ BAT 
or IB I S-ISGRI (sensitive above 15 keV; iBarthehnv et al.l 
120051 : iLebrun et al.1 12003L respectively) are able to de- 
tect only the hottest (fcTbb i^l-5-2 keV) and relatively 
long-lived phases of bursts (see Figure [5]). Thus in 
general the radiated energy and duration measured in 
the hard X-ray band (i? bs and i bs, respectively) are 
lower limits to the total burst fluence and duration (Eh 
and th)- This was illustrated by the long burst from 
XTE J1701-407 detected with both the XRT (sensitive 
in the 0.5-10 keV band) a nd BAT instruments onboard 
Swift (jLinares et al.l l2009h . In that case -Eb and th ex- 
ceeded those measured in the BAT (15-50 keV) band by 
a factor of Eh I E Q hs —3.9 and ibAobs —20, respectively 
(jLinares et al.ll2009[) . 

In order to quantify this combined bandpass and back- 
ground effect for the GBM bursts from 4U 0614+09 and 
to estimate the total radiated energy and duration, we 



simulated GBM light curves using two of the bursts re- 
ported by (KKj): the long PRE burst shown in Figured] 
and a normal duration burst. We created a response ma- 
trix for one GBM Nal detector (nO, 14° off-axis; CTIME 
mode) and used the spectral parameters from both bursts 
(with a 10 s time resoluti on) to simulat e net light curves 
within Xspec (v. 12.7.0: lArnaudl 119961 ). We calculated 
count rates in the energy range 12-25 keV (CTIME chan- 
nel 1), the band used in our systematic search fSec. 12.11) . 
We then added the net burst counts to a real background 
light curve which includes contributions from particle 
background and X-ray sources in the field of view, both 
variable typically on timescales >100 s. Finally, we mea- 
sured the energy (from the integrated fluence, using the 
3.2 kpc distance) and duration of each burst in the sub- 
interval detected by GBM. 

The resulting simulated light curves are shown in Fig- 
ure [6] (right panels). The long burst seen in the GBM 
band would present a characteristic slow rise (i r i S e— 65 s) 
and last for a total of about 140 s. The simulated normal 
duration burst, on the other hand, has a rise time (6.4 s), 
energy (3.6 xlO 39 erg) and duration (30 s) similar to all 
GBM bursts presented herein (Table [2]) . This is shown 
in Figure |5] by comparing two of the GBM bursts (left 
panels; B4 & B8) with the simulated long and normal 
burst light curves (right panels). We conclude that all 
bursts detected from 4U 0614+09 during the first year of 
our GBM X-ray burst search are normal duration bursts, 
yet as discussed below they are more energetic than most 
normal duration bursts. 

Table [3] presents the duration and energy of the long 
and normal duration bursts, as measured with three dif- 
ferent instruments and bandpasses: GBM, FREGATE 
and PCA. The GBM values are estimated from the sim- 
ulations explained ab ove, while FREGATE and PCA 
values are taken from IKIOL The bolometric correction 
factors for energy and duration are given for each case. 
FREGATE was able to detect a slightly longer portion 
of the bursts. Although FREGATE also consisted of 
Nal scintillation detectors with effective area similar to 
that of the GBM detectors, its lower energy threshold 
(6 keV) and smaller field of view (which reduces X-ray 
background) made it sensitive to a larger fraction of the 
bursts. 

The bandpass used has a clear effect on the measured 
burst energy, with Eh/E D h s in the range 1. 3-1.7 (1.3-4 if 
we in clude the XRT/BAT long burst from lLinares et ahl 
2009). Wc note that the same is true for the commonly 
used characteristic burst timescale -Eb/Lpeak- The im- 
pact on the measured burst duration is much more dras- 
tic, as th/t hs can be anywhere between 1.5 and 55. This 



480 



> 

CD 

OJ 
C\J 

in 

CD 

-i—' 

03 
i 

O 

o 

~B 
o 



320 
560 



~ 400 



320 



" - B4(n1) , 

. , , , i , , , , 


■ i i i | i i i i 
.... i ... . 


- B8(n3) = 

aJLsjfeafc "iVfc^ 




. . . . i .... i .... i ... . 



720 
640 
560 
480 
400 
320 
240 
560 

480 

400 

320 

240 



Long K10 



Normal K10 



-100 100 
Time (seconds) 



-100 100 
Time (seconds) 



Fig. 6. — Right panels: Simulated GBM light curves in the 12-25 keV band at 10 s time resolution for a pointing offset of 14°, using the 
spect ral ev olution measured by HETE-FREGATE during a long (top) and normal duration (bottom) burst from 4U 0614+09 (both taken 
from IK 1(31) . Left panels: Light curves for two of the GBM bursts from 4U 0614+09 reported herein, in the 12—25 keV band with 8 s time 
resolution and including background rate. In both cases the pointing offset was 14°; detector number is indicated between parenthesis next 
to the burst ID (Tables [T] and EJ . 

is because cooling ta ils can extend for a long time (e.g.; 
lin't Zand et al.ll2009f) and are detected only with soft X- 
ray detectors, yet most of the energy is radiated dur- 
ing the hottest, most luminous phases of bursts which 
are captured to a greater extent by hard X-ray detec- 
tors. Thus in general burst energy constitutes a more 
robust (less band-dependent) observable than burst du- 
ration. As normal duration bursts from 4U 0614+09 
have not been observed with pointed soft X-ray detec- 
tors, we cannot find the bolometric correction factors for 
the case of normal duration bursts from 4U 0614+09 ob- 
served with GBM. Based on our GBM simulation and 
the FREGATE normal duration burst, we estimate that 
in this case E^/E^s —2 (Table [3]), i.e., normal bursts are 
about twice as energetic as seen in the GBM band. 



4.2. 4U 0614+09 and low-rh bursts 

In Figure [7] we compare the i rcc measured in 
4U 0614+09, when rh was between 0.6% and 1.3%, to 
the constraints available for other low-rh bursters and 
a few selected systems. The t rcc -rh relations for two 
"well-behaved" bursters accreting at a few percent of 
tti-fm are shown fo r comparison (| Galloway et al.l 120041 : 
iFalanga et all 120111 ). In those cases t rec is shorter than 
a day and thus much easier to measure. For compar- 
ison, He ignition curves for pure He accretion are also 
shown, for two different values of the heat flux from 
the crust (Qh=l MeV/nucleon and Qb=3 MeV/nucleon; 
see iCumming fc Bildstenl 120001 for details) . The large 
change in predicted i rcc illustrates the sensitivity of low- 
rh bursts to the thermal properties of the crust. These 
models, however, predict burst energies much higher 
than the observed ones. Detailed modeling of the burst 
properties is beyond the scope of this work, and we refer 



the reader to IKIOI for further discussion. 

It is interesting to compare the distribution of burst en- 
ergies from 4U 0614+09 with that of the normal duration 
and long/intermediate burst populations (see Figure |8]). 
As discussed in Section 14. 1[ burst energy is less sensi- 
tive to bandpass effects than burst duration, and may 
constitute a more robust observable for burst classifica- 
tion. 4U 0614+09 has on average larger E^ than normal 
duration bursts, yet smaller than long/intermediate 
duration burs ts. T he energies of the two long bursts from 
4U 0614+09 (jKlOl ). however, are consistent with the rest 
of the long burst population. Even after bandpass effects 
are corrected for, 4U 0614+09 bridges the gap between 
normal duration and long/intermediate bursts (Fig. [8]). 
This could be due to a fuel composition different than 
most ultracompact LMXBs that would give rise to atyp- 
ical bursts, if as it has been suggested the donor star has 
unusu al abundances (H and He defficient; iWerner et al.l 
2006, see also Sec.QJ. 

Alternatively, the previously observed bimodal distri- 
bution of burst durations could be due to a selection 
effect. Prior to our GBM campaign, observations of low- 
rh bursts were mostly based on serendipitous detections 
with hard X-ray wide-field instruments. If such serendip- 
itous and so far relatively scarce detections are biased 
towards the longest and most energetic events, a popu- 
lation of infrequent, low-rh bursts with energies and du- 
rations between normal and long bursts could have gone 
unnoticed. Our GBM all-sky search (Sec. 12. ip is the first 
systematic survey of low-rh bursts and is less biased to- 
wards long and energetic events. A continuous distribu- 
tion in E\y (between 10 39 erg and 10 41 erg) is actually 
not surprising from the standpo int of He ignition mod- 
els (Cumming Bildsten 2000). This is, as far as we 
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to the shortest wait time between bursts that we find, 2.8 d, and the horizontal error bars sh ow the range in L veIB and rh. Open (gray 
dashed) rectangles show the la confidence intervals obtained from 6 yr of WFC data, taken from in't Zand ct al. (2007, only systems below 
50% -^Edd are shown). The solid o pen black rectangle s hows the WFC limits on t Tec for 4U 0614+0 9. Open squares and triangles show 
tree measurements in GS 1826-24 llGailoway et all 120041 ) and I GR J17511-3057 ijFalanga et al.1120111 '1. respectively. The solid black line 
shows the t rcc cc 1/m relation found by Ga lloway et al.l 1120 041) for G S 1826-24. Green dashed and purple dotted lines show the t rec (m) 
relation predicted by He ignition models (Cummins & Bildstcn 2000, for pure He accretion and solar metalicity) for a base heat flux of 
1 MeV/nucleon and 3 MeV/nucleon, respectively. The arrow below the bottom axis shows the change on the Eddington fraction if a value 
of L E dd=3.8xl0 38 erg s" 1 is used (Sec. [21}. 



know, the first time that such continuous distribution is 
observed. Extending this study to similar bursters will 
determine if the distribution of burst energies that we 
find (Figure [5J is unique to 4U 0614+09 or a common 
feature of low-m bursters. 

The average burst energy in our sample, after the bolo- 
metric correction factor of ~2 is applied (Sec. I4.1j) . is 
5xl0 39 erg. For the measured < rcc (12 d) and L pcrs 
(2.3xl0 36 erg s _1 ), this corresponds to an accretion-to- 
burst fluence ratio a ~ 560, while the accreted column 
depth at the average inferred rh (1.1 xlO 3 g cm~ 2 s _1 ) is 



y ~ 1x10 g cm" , similar to the values obtained from 
light curve fitting (2xl0 9 g cm~ 2 for the normal dura- 
tion burst in lKlOl ). However, complete ignition of such a 
thick layer would give rise to more energetic bursts than 
observed, assuming the ener gy re lease expected for pure 
He burning (as discussed bv lKlOl ). For the shortest f wait 
burst (B13), using the average L pcis and rh during Jan- 
uary 2011 (1.8xl0 36 erg s" 1 and 0.9xl0 3 g cur 2 s" 1 , 
respectively) yields a similarly high value of a ~506. This 
indicates that even the shortest t rcc bursts in 4U 0614+09 
are sub-energetic. Thus after our improved t lec measure- 
ments and considering bolometric/bandpass corrections 



to -Eb, the question of why are most bursts in 4U 0614+09 
sub-energetic remains open. We also note that the short- 
est i wa it burst (B13) has the lowest E\> in our sample, as 
expected if the mass accumulated before ignition is also 
lowest. We found the period with the shortest i rcc near 
the minimum in rh (Sec. [3]). This behavior is opposite 
to what ignition models predict and it might be due to 
different rates of stable H (CNO) burning operating at 
slightly different rh (although the apparent lack of H in 
the fuel makes this scenario unlikely). Alternatively, it 
could reflect changes in the L peTS -rh relation (i.e., in the 
radiative efficiency of the accretion flow) unaccounted 
for. Continued monitoring may determine whether this 
is a weeks-long fluctuation in the observed burst rate or 
it constitutes a systematic behavior of 4U 0614+09. 

As can be seen in Figure [71 our GBM campaign in its 
first year has improved the accuracy of i rec measurements 
at low rh, placing new constraints on models of low-m 
ignition. Most estimates available to date were based 
on BeppoSAX-WFC observations and provided mostly 
order-of-magnitude constraints on t rcc (Fig. [Jj. Extend- 
ing our GBM measurements to other low-rri bursters and 
comparing them to theoretical models, considering po- 
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Fig. 8. — Bottom: Burst peak luminosity vs. energy, for the GBM bursts from 4U 0614+09 presented herein (filled red squares; error bars 
show estimated bolometric correction factor, see Sec. I4.lt . for all kno wn bursts from 4U 614+09 including those reported in KlTj (open 
red squares), for all bursts f rom the RXTE cata log (filled gray circles, Galloway et al. 20081) . for the 11 Hz pulsar in the globular cluster 
Terzan 5 (open gray circles, Linares et al. 2012) and for the long/intermediate bursts from the MINBAR catalog (filled blue triangles). 
Middle and top panels: Histograms showing the number of bursts per energy bin for the diffe rent known samples, as indicated. 4U 0614+09 
bridges the gap between the "normal" and "long/intermediate" burst populations fSec 14.2b . 



tential changes in fuel composition across bursters, can 
make precise measurements of crustal heat flux possible. 

5. SUMMARY 

• We have presented 15 normal duration thermonuclear 
bursts from 4U 0614+09 observed with Fermi-GBM be- 
tween March 2010 and March 2011, when the source 
was accreting near 1% of the Eddington rate. We mea- 
sured a burst recurrence time of 12±3 d and find the 
closest burst pair seen from 4U 0614+09 to date, 2.8 d 
apart. 

• We find no long/intermediate bursts in one year of 
data, and set a lower limit to their recurrence time of 
62 d. The 2001-2002 "calm" (low long-term variabil- 
ity) period in the persistent soft X-ray emission during 
which long/intermediate bursts were detected has not 
recurred. 

• We quantify the bandpass effect on the observed burst 
duration and energy. We find that the total burst du- 
ration can be 1.5-55 times longer than that measured 



in the hard X-ray band, whereas this bolometric correc- 
tion factor is 1.3-1.7 for burst energy (which is therefore 
a more robust, band-independent observable). 

• The burst energies in 4U 0614+09, between 8xl0 38 erg 
and lxlO 40 erg after our bolometric correction is ap- 
plied, overlap between the energies of the normal and 
long burst populations. This continuous distribution 
in burst energy provides new observational evidence 
that long/intermediate bursts are an extreme case 
of low- m bursts. We suggest that the apparent bi- 
modal distribution of durations that defined normal and 
long/intermediate bursts could be due to a selection ef- 
fect if only the longest and most energetic low-m bursts 
could be detected to date. 
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